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ABSTRACT 

Selecting centrally quiescent galaxies from the Sloan Digital Sky Survey (SDSS) to create high signal-to- 
noise (> 100 A -1 ) stacked spectra with minimal emission line contamination, we accurately and precisely 
model the central stellar populations of barred and unbarred quiescent disk galaxies. By splitting our sample 
by redshift, we can use the fixed size of the SDSS fiber to model the stellar populations at different radii within 
galaxies. At 0.02 < z < 0.04, the SDSS fiber radius corresponds to w 1 kpc, which is the typical half-light radii 
of both classical bulges and disky pseudobulges. Assuming that the SDSS fiber primarily covers the bulges at 
these redshifts, our analysis shows that there are no significant differences in the stellar populations, i.e., stellar 
age, [Fe/H], [Mg/Fe], and [N/Fe], of the bulges of barred vs. unbarred quiescent disk galaxies. Modeling the 
stellar populations at different redshift intervals from z = 0.020 to z = 0.085 at fixed stellar masses produces an 
estimate of the stellar population gradients out to about half the typical effective radius of our sample, assuming 
null evolution over this ss 1 Gyr epoch. We find that there are no noticeable differences in the slopes of the 
azimuthally averaged gradients of barred vs. unbarred quiescent disk galaxies. These results suggest that bars 
are not a strong influence on the chemical evolution of quiescent disk galaxies. 

Subject headings: galaxies: abundances — galaxies: stellar content — galaxies: structure — galaxies: evolu¬ 
tion 


1. INTRODUCTION 

The stellar populations of galaxies offer a unique view of 
the formation and evolution of galaxies. In this work, we use 
stellar populations to explore the evolution of barred quiescent 
disk galaxies. 

Bars are predicted to affect the properties of their host 
galaxies through the slow rearrangement of energy and mass, 
a pro c ess called secular evo l ution dKormendv & Kennicutti 
120041 lAthanassoulal 1201 3bt ISellwoodl 120141) . This - re¬ 
distribution of material leads to many effects, including 
the flattening of the i nitial global abun d ance gradients 
of th e ir host galaxies (IFriedli et all 11 994t iFriedli & Benzl 


19951 iMartinet :& FriedliT *1997|: Minchev & Famaevl 120101 


Di Matteo et al.1 [2013tlKubryk et al.lT2013lk which results in 
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bulge abundances that are dependent on the presence of 
star formation. Specifically, if star formation is present, 
then the gaseous metallicities of bulges in barred galax¬ 
ies are predicted to increase compared to unbarred galaxie s 
(IFriedli et al.ll 1991 IFriedli & Benll995l: iMartel et alJflOUl) . 

while the bulge stellar metallicities stay relatively unchanged 
(IFriedli et al.ll l994b . In this case, bars are also predicted to 
produce younger stellar populations at the center and the ends 
of the bar structure (IWozniakl 120071) . But if star formation 
is absent, then both the gaseous and stellar metallicities of 
bulges in barred galaxies are predicted to decrease due to di¬ 
lution by lower metallicity stars and gas from further galacto- 
centri c distances, assuming an initial negative metallicity gra¬ 
dient (IFriedli et alJl!994llDi Matteo et alJl2013l) . 

So far, works that have studied the bulge gas-phase metal¬ 
licities of barred galaxies have only considered star-forming 
galaxies, which means that the bulge gaseous metallicities 
of barred galaxies are expected to be enhanced relative to 
unbarred galaxies. While there are works that agree with 
this prediction, there are others that do not. For example, 
lEllison et al.1 ( 201 D. who used a large sample of galaxies 
from the Sloan Digital Sky Survey (SDSS), found that the gas- 
phase metallicities in the inner few kpc of barred galaxies are 
higher (by ~0.06 dex) than unbarred galaxies at a given stellar 
mass. But ICacho et all (120141) . who also used a large sample 
of SDSS galaxies, found that for a given stellar mass there 
are no significant differences in the gas-phase metallicities at 
the inner few kpc of ba rred and unbarred galaxies (see also 
iHenrv & Wo rlhev 1 999!) . However, in contrast to all these re¬ 
sults. iDutil & Ro\T i 1999) and IConsidere et all (120001) found 
that the central gas-phase metallicities of barred galaxies are 
lower than that of unbarred galaxies. 

Works that have studied the bulge stellar populations of 
barred galaxies have considered both star-forming and qui¬ 
escent galaxies. However, in analyzing their results, most 
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of these works did not separate star-forming galaxies from 
quiescent galaxies, thus producing results that are difficult 
to interpret. Nonetheless^ we summarize their findings for 
the benefit of the reader: iMoorthv & Holtzmanl (i2006l) found 
that barred galaxies have higher bulge stellar metallicities 
than unbarred galaxies at a given velocity dispersion, which 
is consistent with the findings of lPerez & Sanchez-Blazauezi 
(1201 lh . who additionally found that the bulges of a sample 
of barred SO galaxies are slightly a-enhanccd compared to 
that of unbarred SO galaxies. These results, however, are in 
contrast to t hose oflC oelho & Gadottil (1201 11) . I Williams et alJ 
(120121) . and lCacho etal.1 (120141) : they all found no differ¬ 
ence in the stellar metallicity of bulges. Moreover, there 
is a similar disagreement about the bulge_stellar_ages of 
barre d and unbarred galax ies: IPerez & Sanchez-Blazauez 
1201 1L iWilliams et alj[2 0l2l and iC acho et alJ 120141 found no 
difference, while iCoelho & Gadottil (1201 11) did. 

Clearly, there is no consensus on whether bars affect the 
bulge abundances of their host galaxies. However, the lack of 
any strong bulge abundance differences in these works sug¬ 
gests that any potential effects from bars are probably small. 

Turning to global gradients, which are predicted to be flatter 
in barred galaxies regardless of the presence of star formation, 
reveals a similarly confused landscape. While some studies 
found that the gas-phase metallicity gradients of barred galax- 


1992; Martin & Roy 

1994| Zaritsky et al. 

1994), other studies 

found no difference 

Sanchez et al. 2014| 

Ho et alj2()15). 


Studies on the stellar population gradients of barred galax¬ 
ies have also found conflicting results^ Using slits placed 
along the bars of 20 galaxies, IPerez et ahl d2009l) found a 
variety of stellar metallicity and age gradients, including 
positi ve, negative, and null gradients_ (see also Perez et al 
12007 ). For two of these galaxies, ISanchez-Blazquez et al 


( 201 1) found that the stellar metallicity and age gradients 
are flatter along the bar than along the disk. In agreement, 
IWilliams et al.l 00121) used a sample of 28 galaxies to show 
that the stellar metallicity gradients of boxy/peanut-shaped 
bulges (which are presu mab ly, indicative of the presence 
of a bar; lAthanassoulal [2005 ) are shallower than that of a 
sample of unbarred early-type galaxies. But most recently, 
using 62 face-on sp ir al galaxies from the CALIFA survey 
(ISanchez et al.l 120121) . ISanchez-Blazquez et alJ d2014l) found 
no difference in the stellar metallicity gradient or age gradient 
(their gradients extend out to the disk region) between barred 
and unbarred galaxies. 

Similar to bulge abundances, there is no consensus on 
whether bars affect the abundance gradients of their host 
galaxies, indicating that any potential effects due to bars are 
small. Therefore highly accurate and precise chemical analy¬ 
sis is needed to robustly detect them. 

This work aims to achieve this goal by using very high 


signal-to-noise (S/N >100 A -1 ) stacks of barred and unbarred 
galaxies, identified by Galaxy Zoo 2 dWillett et al. 112013l) . to 
model their stellar populations. Aside from the unprecedented 
S/N stacks of barred galaxies, an additional important and 
unique feature of our study is the minimization of emission 
lines by selecting only quiescent galaxies. This latter feature 
circumvents the need to model the emission line contribution, 
which is complex and introduces more uncertainty to the re¬ 
sulting stellar population parameters dConrovl 120131) . More¬ 
over, with quiescent galaxies, we are comparing to predictions 
that have not been investigated, i.e., simulations without star 


formation. 

We describe the data, sample selection, stacking, and stel¬ 
lar population synthesis model in §2 and §3. §4 presents our 
results, which indicates that there are no significant differ¬ 
ences in the stellar populations of barred and unbarred qui¬ 
escent disk galaxies. We present a comparison to past works 
in §5 and discuss our results in §6. Finally, we close in §7. 
Throughout this paper, we assume a flat cosmological model 
with //(> = 70 km s -1 Mpc -1 , Q„, = 0.30, and Ua = 0.70, and 
all magnitudes are given in the AB magnitude system. 

2. DATA 

Since our study aims to model the stellar populations of 
barred and unbarred quiescent disk galaxies, we must procure 
spectra of a large sample of these types of galaxies. Thus we 
use the spectroscopic Main Galaxy Sample in the Leg acy area 
of the SDSS Data Release Seven (SDSS DR7; iStrauss et alJ 
l2002tlAbazaiian et alj|2009l) : these spectra are measured with 
a 3" fiber, hence the results presented in this work only con¬ 
cern the central regions of our galaxies. 

2.1. Sample selection 

In order to most accurately and most precisely model the 
detailed abundance patterns of galaxies, passive galaxies with 
little-to-no emission lines are required. Thus we select quies¬ 
cent galaxies as described by iPeek & Gravesl (120101) . i.e., no 
Ha or [O II ] A3727 emission; these criteria have been shown 
to exclude star-forming galaxies and AGN hosts ([Graves et all 
120091) . Compared to a red sequence defined by specific star 
formation rate or rest-frame g—r, this spectroscopic quiescent 
criteria is much more strict—about 50% of the galaxies in the 
red sequence make it into our quiescent sample. However, 
the fraction of face-on disk galaxies and the fraction of barred 
galaxies, which are defined below, are similar in both the red 
sequence and our spectroscopically-defined quiescent sample, 
indicating that our sample selection is not picking a special 
subset of red sequence galaxies. The presence of quiescent 
disks is consistent with recent work that showed that galax¬ 
ies with disks have made up a significant fraction of the re d 
sequence since z ~ 1 dBundv et all2 010; Ma sters et al.ll2010l) . 

In addition to our quiescent selection, we follow the SDSS 
DR7 caveanFl and select galaxies with median S/N>10 A -1 
and with velocity dispersions of 70 km s -1 < a < 420 km s -1 . 
For these selections, we have utilized the MPA-JHU DR7 
value-added catalogsQ 

2.1.1. Selecting barred and unbarred disk galaxies 

To select barred disk galaxies from our quiescent sample, 
we first selec t face-o n di sk galaxies ba s ed on the criteria of 
iMasters et al.l (1201 If) and ICheung et alJ (120131) . Namely, we 
select galaxies with axis ratios (b/a) larger than 0.5 since it 
is difficult to identify bars in edge-on galaxies (b/a measure¬ 
ments are from the GIM2D single Sersic fits bv lSimard et all 
1201 ll) . We then use the debiasecQ Gal axy Zoo 2 (GZ2) vi- 
sual morphological classifications from IWillett et al.l d201 3l) 

13 http://classic.sdss.org/dr7/products/spectra/index.html#caveats 

14 http://www.mpa-garching.mpg.de/SDSS/DR7/ 

15 The debiasing procedure accounts for the deterioration of the image 
quality due to increasing distance of galaxies, i.e., with increasing redshift. 
The underlying assumption is that galaxies of a similar luminosity and size 
will share the same average mix of morphologi es w ithin the epoch covered 
by our redshift range (0.020 < z < 0.085; see IWillett et al.1120131 for more 
details). 
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Figure 1. Graphic summary of this study. Stellar mass (M*) vs. redshift (z) of our quiescent galaxy sample is shown in the left-center panel. Light-blue hatches 
mark the galaxies used to analyze the bulge stellar populations, and the light-green hatches mark the galaxies used for the stellar population gradient analysis. 
The spectra of the barred and unbarred galaxies in the black rectangle at 0.02 < z < 0.04 and at 10.25 < log M* /Mq < 10.40 are shown at the top and bottom, 
respectively. Each individual spectrum is plotted in grey, and the stacked spectra of the barred and unbarred galaxies are plotted in their respective colors. These 
stacks have the high signal-to-noise (S/N > 100 A~*) needed to obtain accurate and precise stellar population parameters. Example images of barred and unbarred 
galaxies from the highlighted bin are also displayed, along with their DR7 SDSS object IDs. Finally, the center-right panel shows the ratio of the barred and 
unbarred stacked spectra at several Lick indices. This spectral comparison illustrates the extremely subtle differences that are analyzed in this work, and hence 
the need for high S/N and full-spectrum modeling. 


to select non-edge-on disk galaxies by requiring that for each 
galaxy, at least a quarter of its classifications involved answer¬ 
ing the bar question. 

A relevant aside is to compare our quiescent disk sample to 
the red spiral sample of lMasters et al.lfl2010t) : we find that less 
than 1% of their sample make it into our sample. The main 
reason is our requirement of little-to-no Ha or [O II] A3727 
emission, which is app arently present in almost all red spiral 
galaxies (iMasters et alj d2010l) also excluded large bulges and 
imposed a more strict inclination cut). This comparison indi¬ 
cates that the disks of this sample of galaxies are distinctive in 
that they have little to no detectable spiral structure, i.e., they 
are featureless disks. 

From this face-on quiescent disk sample, we select barred 
galaxies with a debiased bar vote fraction threshold of 0.5, 
i.e., galaxies with pbm > 0.5 are considered barred. This 
bar threshold has been shown to be a reliable indicator of 
strong bar features—almost all galaxies with /; h ar > 0.5 were 
classifie d as possessing a strong ba r by [N air & A brahaml 
(1201 Obi) (see IMasters et alJ IMJ and IWillett et al.l l2013h — 
and has been adop t ed by several past Galaxy Zoo works 
(Mast ers et al.11201 it IMasters et al.H2012t iMelvin et al.l 120141 : 

iCheung et al.l 12015b . Using a more strict debiased bar vote 

fraction threshold of pb ar = 0.8 does not change our qualitative 
results. We present representative images of barred galaxies 
in a given stellar mass and redshift interval in Fig. 0 to illus¬ 
trate the reliability of the GZ2 visual classifications. 


To select unbarred disk galaxies, we choose galaxies with 
a debiased bar vote fraction of zero, i.e., p\ YM = 0.0, from 
the face-on quiescent disk sample described above. This 
bar vote fraction choice is different from previous Galaxy 
Zoo works because we want to minimize contamination from 
weak bars, which h as been shown by M asters et all (1201 2l) and 
Will ett et al.l (120131) to likely have 0.1 < pw M - < 0.5. 

The selection of face-on unbarred quiescent disks is diffi¬ 
cult and there is no perfect technique. However, comparing 
the structural parameters (global Sersic index, central surface 
stellar mass density, central velocity dispersion, and bulge-to- 
total ratio) of our barred and unbarred quiescent disk samples 
for a given stellar mass bin shows a similar distribution, indi¬ 
cating that we select similar types of galaxies (see Figs. [6]and 
[7]in Anncndix[/\li. Moreover, comparing the structural param¬ 
eters of our quiescent disk samples to a sample of quiescent 
spheroid-dominated galaxies shows differences at the lowest 
stellar masses, indicating that our samples are not severely 
contaminated by spheroid-dominated objects. 

And finally, since our study concerns isolated galaxies, 
we eliminate merging galaxies by discarding all galaxies 
with a_ Galaxy Zoo merging parameter, p mg , larger than 0.4 
dDarg et al.|l2010i) . 

2.1.2. Quiescent barred galaxies 

Although barred galaxies are commonly thought of as star¬ 
forming, recent works have shown that barred galaxies are 
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Table 1 

Properties of Stacked Spectra 



z Bin 

M* Bin 

Number” 

S/N b 



(Mq) 


(A' 1 ) 

Bulge-barred 1 

0.020 < z < 0.040 

9.60 < log M „ < 9.90 

13 

80 

Bulge-barred 2 

0.020 < z < 0.040 

9.90 < log M» < 10.10 

32 

130 

Bulge-barred 3 

0.020 < z < 0.040 

10.10 < log M» < 10.25 

39 

170 

Bulge-barred 4 

0.020 < z < 0.040 

10.25 < log M, < 10.40 

62 

257 

Bulge-barred 5 

0.020 < z < 0.040 

10.40 < log M, < 10.55 

71 

307 

Bulge-barred 6 

0.020 < z < 0.040 

10.55 < log M, < 10.70 

66 

307 

Bulge-barred 7 

0.020 < z < 0.040 

10.70 < log M, < 10.85 

26 

189 

Bulge-barred 8 

0.020 < z < 0.040 

10.85 < log M, < 11.10 

10 

142 

Bulge-unbarred 1 

0.020 < z < 0.040 

9.60 < log M* < 9.90 

39 

133 

Bulge-unbarred 2 

0.020 < z < 0.040 

9.90 < log M, < 10.10 

81 

216 

Bulge-unbarred 3 

0.020 < z < 0.040 

10.10 < log M* < 10.25 

84 

253 

Bulge-unbarred 4 

0.020 < z < 0.040 

10.25 < log M* < 10.40 

83 

294 

Bulge-unbarred 5 

0.020 < z < 0.040 

10.40 < log M* < 10.55 

87 

317 

Bulge-unbarred 6 

0.020 < z < 0.040 

10.55 < log M, < 10.70 

40 

248 

Bulge-unbarred 7 

0.020 < z < 0.040 

10.70 < log M, < 10.85 

18 

178 

Bulge-unbarred 8 

0.020 < z < 0.040 

10.85 < log M* < 11.10 

10 

141 


Gradient-barred 1 

0.020 

< z 

< 

0.040 

10.10 

< 

log 

M„ 

< 

10.40 

101 

306 

Gradient-barred 2 

0.040 

< z 

< 

0.055 

10.10 

< 

log 

M« 

< 

10.40 

131 

272 

Gradient-barred 3 

0.055 

< z 

< 

0.067 

10.10 

< 

log 

M, 

< 

10.40 

59 

153 

Gradient-barred 4 

0.020 

< z 

< 

0.040 

10.40 

< 

log 

M, 

< 

10.70 

137 

435 

Gradient-barred 5 

0.040 

< z 

< 

0.055 

10.40 

< 

log 

M* 

< 

10.70 

139 

338 

Gradient-barred 6 

0.055 

< z 

< 

0.067 

10.40 

< 

log 

M* 

< 

10.70 

143 

261 

Gradient-barred 7 

0.067 

< z 

< 

0.077 

10.40 

< 

log 

M* 

< 

10.70 

148 

235 

Gradient-barred 8 

0.077 

< z 

< 

0.085 

10.40 

< 

log 

M, 

< 

10.70 

66 

158 

Gradient-barred 9 

0.020 

< z 

< 

0.040 

10.70 

< 

log 

M, 

< 

11.00 

36 

233 

Gradient-barred 10 

0.040 

< z 

< 

0.055 

10.70 

< 

log 

M t 

< 

11.00 

90 

318 

Gradient-barred 11 

0.055 

< z 

< 

0.067 

10.70 

< 

log 

M * 

< 

11.00 

92 

262 

Gradient-barred 12 

0.067 

< z 

< 

0.077 

10.70 

< 

log 


< 

11.00 

86 

213 

Gradient-barred 13 

0.077 

< z 

< 

0.085 

10.70 

< 

log 

Mt 

< 

11.00 

94 

207 

Gradient-unbarred 1 

0.020 

< z 

< 

0.040 

10.10 

< 

log 

Mt 

< 

10.40 

167 

383 

Gradient-unbarred 2 

0.040 

< z 

< 

0.055 

10.10 

< 

log 

Mt 

< 

10.40 

509 

527 

Gradient-unbarred 3 

0.055 

< z 

< 

0.067 

10.10 

< 

log 

Mt 

< 

10.40 

341 

373 

Gradient-unbarred 4 

0.020 

< z 

< 

0.040 

10.40 

< 

log 

Mt 

< 

10.70 

127 

399 

Gradient-unbarred 5 

0.040 

< z 

< 

0.055 

10.40 

< 

log 

Mt 

< 

10.70 

429 

604 

Gradient-unbarred 6 

0.055 

< z 

< 

0.067 

10.40 

< 

log 

Mt 

< 

10.70 

687 

623 

Gradient-unbarred 7 

0.067 

< z 

< 

0.077 

10.40 

< 

log 

Mt 

< 

10.70 

842 

600 

Gradient-unbarred 8 

0.077 

< z 

< 

0.085 

10.40 

< 

log 

Mt 

< 

10.70 

441 

411 

Gradient-unbarred 9 

0.020 

< z 

< 

0.040 

10.70 

< 

log 

Mt 

< 

11.00 

26 

213 

Gradient-unbarred 10 

0.040 

< z 

< 

0.055 

10.70 

< 

log 

Mt 

< 

11.00 

152 

418 

Gradient-unbarred 11 

0.055 

< z 

< 

0.067 

10.70 

< 

log 

Mt 

< 

11.00 

292 

493 

Gradient-unbarred 12 

0.067 

< z 

< 

0.077 

10.70 

< 

log 

Mt 

< 

11.00 

423 

509 

Gradient-unbarred 13 

0.077 

< z 

< 

0.085 

10.70 

< 

log 

Mt 

< 

11.00 

472 

484 


a Total number of galaxies in the stacked spectrum, 
b Effective median S/N of the stacked spectrum. 


actually m ore frequent among the red sequence disks than the 
blue c loud dMasters et a 1 .11 201 lb iLee et alJ2012UCheung et all 
120131) . And when con sidering absolute nu mbers in the 
volume-limited sample of lCheung et all (120131) . the numbers 
of red and blue barred galaxies are comparable. Finally, al¬ 
though our spectroscopically-defined quiescent sample is ap¬ 
proximately half of the photometrically-defined red sequence, 
the bar fraction of the two samples is similar, indicating that 
our quiescent barred galaxy sample is not an unusual class of 
red sequence barred galaxies. 

2.2. Stacking 

Because detailed abundance analysis requires high S/N 
spectrum (S/N > 100 A -1 ; ICardiel et al.lU9^8l) . we need to 
stack many SDSS spectra (each SDSS spectra has a typical 


S/Nw 20 A ') to achieve the necessary S/N (e.g. JGraves et all 

l200l . 

We use two stacking schemes that correspond to our two 
science goals: (1) bins of stellar mass at 0.02 < z < 0.04 for 
our bulge analysis and (2) bins of redshift at three mass in¬ 
tervals for our gradient analysis (see Table [I] galaxy stellar 
masses [not fiber stellar masses] are taken from the MPA-JHU 
DR7 value-added catalogs); we explain the motivations for 
these schemes in ^4] 

Each individual spectrum was continuum-normalized using 
gaussian convolution, smoothed to an effective velocity dis¬ 
persion of 300 km s -1 , corrected for gala ctic extinction us¬ 
ing the parameterization of ICardelli et alj (119891) . and flux- 
calibrated a ccording to I Yanl (1201 lh . The empirical flux cali¬ 
bration from lYan ( 2011 ) results in relative spectrophotometric 
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precision of 0.1%, increasing the accuracy and precision of 
our model fitting, which we describe in the next section. Each 
spectrum contributes equally to a stack, but portions of spec¬ 
tra that are 5 ct deviations from the initial stack are trimmed. 
Only a small percentage of spectra are trimmed in this man¬ 
ner, and the exclusion of this trimming does not affect our 
final results. The typical S/N of our stacks is > 100 A -1 ; see 
TableQ]for each stack’s S/N. 

We compare the stacks of our barred and unbarred samples 
in Fig. [8] and [9] of Appendix [B] Focusing on the ratio of the 
stacks at the Fick indices, we find that there are no signifi¬ 
cant differences, which is supported by Fig. [10] which plots 
several Fick indices as a function of mass for barred and un¬ 
barred quiescent disk galaxies. This initial analysis supports 
our findings using full-spectrum fitting in $4] 

The total number of barred and unbarred galaxies that make 
up the stacks for our bulge analysis is 319 and 442, respec¬ 
tively. The total number of barred and unbarred galaxies that 
make up the stacks for our gradient analysis is 1,322 and 
4,908, respectively. 


3. MODEL FITTING 


We fit our high S/N stacked spectra with the latest 
version of the stellar population synthesis (SPS) model 
of IConro v & van DokkunJ (1201 2al). which is detailed in 
IConrov et al.1 (12014) (see also Choi et al.l 120141) . We follow 
the methodology of Conrov et al.l (120141) in this present work. 

Briefly, this SPS model uses three sets of stellar isochrones 
and the MIFES dSa nchezzBlazguezetakl 120061) and IRTF 
(ICushing et al.l 12005b iRavner et akT 20091) empirical stellar 
spectral libraries to fit for the full inputted spectrum at four 


wavelength intervals: 4000 A^-700 A, 4700 A-5700 A, 5700 
A-6400 A, and 8000 A-8800 A. The wavelength range 6400 
A-8000 A is not included because it is at the edges of the 
MIFES and IRTF wavelength coverage. A two-part Kroupa 
IMF is assumed. Non-solar abundance patterns are modeled 
with response functions using the ATFAS12/SYNTHE code 
suites (iKuruczl 19701119931) and applied differentially onto the 
template model galaxy spectrum. A Markov Chain Monte 
Carlo (MCMC) algorithm was used to explore the 40 free pa¬ 
rameters of the model, which include: redshift, velocity dis¬ 
persion, two population ages (the age of the dominant popula¬ 
tion and the age of the younger population J3, the mass frac¬ 
tion of the younger population, four nuisance parameterfQ 13 
emission line strengths, the velocity broadening of the emis¬ 
sion lines, and the abundances of C, N, Na, Mg, Si, Ca, Ti, V, 
Cr, Mn, Fe, Co, Ni, Sr, and Ba, and O, Ne, S (the latter three 
are varied in lock-step). The systematic uncertainties (e.g., 
incomplete stellar spectral libraries, unknown aspects of stel- 
lar evolution, metal l icity evolution) a re probably ~ 0.05 dex 
(IConrov et alJl2014t IChoi et afll2014l) . 

We illustrate the quality of our model fits in Figs, fill andfl2l 
of Appendix |Cj it shows that the fits are excellent and more¬ 
over, the quality of the barred and unbarred fits are indistin¬ 
guishable. 


16 Although a quiescent disk galaxy may have a more complex star forma¬ 
tion history, the old ages we derive in ^indicates that this issue should not 
affect our results <Serra & Tragerl20()7l) . 

17 These parameters are meant to describe uncertain aspects of SPS mod¬ 
eling. Two of these are the temperature and the fractional flux from young 
stars or hot horizontal branch stars. One allows for the addition of arbitrary 
amo unts of M giant light, and the final one is a shift in effective temperature 
fsee lConrov & van DokkLiml201 2)4 and lConrov et alJ2014l for more details). 


4. RESULTS 

4.1. Stellar populations of the bulges of barred and unbarred 

galaxies 

To model the stellar populations of galaxy bulges, our anal¬ 
ysis only considers redshifts between 0.02 < z < 0.04 since 
the SDSS fiber radius (l."5) at z = 0.04 corresponds to « 1 
kpc, which is the typical half-light radii of both classical 
bulges and disky pseudobulges dFisher & Drorvl 120101) . The 
lower limit of z = 0.02 corresponds to the redshift where the 
SDSS spectrograph starts to cover [O ll]A3727. The sample 
for our bulge stellar population analysis is highlighted by the 
blue hatches in the log M* — z distribution of our quiescent 
sample in Fig. Q] The stellar masses we consider range from 
log M*/Mg = 9.6 to log M*/M q = 11.1, which we split into 
eight log M* bins; these bins are listed in the top rows of 
Tableland are labeled “Bulge-barred” or “Bulge-unbarred”. 
The barred and unbarred stacked spectra of the highlighted 
bin (0.02 < z < 0.04 & 10.25 < log M»/M 0 < 10.40) are 
shown in the upper and lower panels of Fig. Q] with each 
spectrum plotted in grey. Two images of barred and unbarred 
galaxies from the highlighted bin are also shown. The ra¬ 
tio of the barred and unbarred stacked spectra at several Fick 
indices are highlighted in the center-right panel. There are 
extremely subtle differences in these Fick indices, indicating 
that: (1) high S/N is needed to detect these differences and 
(2) full-spectrum modeling can help bring out features that 
are not covered by Fick indices. Appendix[B]presents a more 
detailed comparison of the stacks and their residuals—they 
show similar features to the comparison in Fig. [7] 

With these high S/N stacks, we model the bulge stellar pop¬ 
ulations of our barred and unbarred quiescent disk galaxies at 
each stellar mass bin. Our results are shown in Fig. [2} it dis¬ 
plays the stellar age, [Fe/H] (a metallicity indicator), [Mg/Fe] 
(a star formation timescale indicator), and [N/Fe] (another star 
formation timescale indicator) vs. stellar mass of barred qui¬ 
escent disks in magenta and unbarred quiescent disks in pur¬ 
ple (this is the color scheme throughout the paper). There is a 
general correlation between all stellar population parameters 
and stellar mass for both populations, similar to the trends be¬ 
tween stellar population parameters and central velocity dis¬ 
persion of early-type galaxies (e.g. JThomas et alll2005i) . 

The errors estimates shown in Fig. [2] are based on the 
full posterior distributions of each fitted parameter from the 
MCMC spectrum -fitting algorithm. Systematic uncert ainties 
are ~ 0.05 dex (IConrov et al.l 120141: IChoi et all 1201 4ft and 
should be assumed in addition to the displayed error esti¬ 
mates. Thus the main result of Fig. [2] is that there are no 
significant differences in the stellar populations of the bulges 
of barred vs. unbarred quiescent disk galaxies, which is sup¬ 
ported by the P-values of > 0.14 (based on the Z-score). 

4.2. Stellar population gradients 

Since the SDSS fiber size corresponds to a larger physi¬ 
cal size at higher redshift, we can probe stellar populations at 
larger physical galactocentric radii by analyzing higher red¬ 
shift galaxies (e.g.. lYan & Blanton! 120121) . Combining this 
approach with the assumption that galaxy properties (e.g., 
galaxy size, bar length, and their abundance patterns) do not 
evolve significantly at a fixed stellar mass over our redshift 
range of 0.020 < z < 0.085, which corresponds to « 1 Gyr, 
we can estimate stellar population gradients. Our upper limit 
of z = 0.085 corresponds to the redshift limit of the debias¬ 
ing procedure applied to the GZ2 vote fractions (IWillett et al.l 
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Figure 2. Bulge stellar populations vs. stellar mass for barred and unbarred quiescent disk galaxies. Namely, a): Age vs. log M», b): [Fe/H] vs. log M», c): 
[Mg/Fe] vs. log A/*, and dj: [N/Fe] vs. log M* at 0.02 < z < 0.04. The error estimates are based on the full posterior distributions of each parameter from 
the MCMC spectrum-fitting algorithm. We display a 0.05 dex floating error bar in each panel as a conservative estimate of the systematic errors, which should 
be assumed in addition to the displayed error bars. There are no significant differences in age, [Fe/H], [Mg/Fe], or [N/Fe] of the bulges of barred and unbarred 
quiescent disk galaxies. 
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Figure 3. Stellar population gradients of [Fe/H], [Mg/Fe], and [N/Fe] at three different stellar mass bins—10.1 < log M*//Wg < 10.4, 10.4 < log M,/Mg < 
10.7, and 10.7 < log M* /Mq < 11.0—for barred and unbarred quiescent disk galaxies. The fourth row superimposes all mass bins. The corresponding redshifts 
are displayed on the top axis. We display a 0.05 dex floating error bar in every panel as a conservative estimate of the systematic errors. There are no clear 
differences in the slopes of gradients of barred and unbarred galaxies. 


120131) . Therefore, we separate our sample into seven bins of 
redshift, from z = 0.020 to z = 0.085, at three stellar mass bins: 
10.1 < log M,/M 0 < 10.4, 10.4 < log M*/Mq < 10.7, and 
10.7 < log M,/M q < 11.0; these bins are listed in the bot¬ 
tom rows of Table Q] and are labeled “Gradient-barred” or 
“Gradient-unbarred.” This sample is highlighted with green 
hatches in the log M* vs. z panel of Fig. Q] 

4.2.1. SDSS fiber size vs. typical galaxy size 

Before we present the stellar population gradients, an im¬ 
portant question that needs to be addressed is “how much of 
our galaxies is covered by the SDSS fiber?” The answer is 


shown in Fig. [4] which plots semi-major half-light radius (r e ; 
from the GIM2D single Sersic fits bv ISimard et a I JlTO 1 111 vs. 
redshift (z). Each pane l dis plays the r e — z distribution of our 
quiescent sample (see 32. lb in the grey contours. Overplotted 
in panels a-c are the r e —z distributions of the three mass bins 
considered in the stellar population gradient analysis. The 
barred and unbarred samples are displayed in solid and dashed 
contours, respectively; there is no evident difference in their 
r e — z distribution. Moreover, there is no clear r e evolution at 
a given mass bin, supporting our assertion of null r e evolution 
in our redshift range. 

Within panels a-c, the colored hatched lines represent the 
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Figure 4. The semi-major half-light radius (r e ) vs. redshift (z) distribution of our quiescent sample is displayed in grey contours in each panel. Panels a-c 
overplots the r e —z distribution of the three mass bins of barred and unbarred galaxies used in the stellar population gradient analysis The radius range that the 
SDSS fiber can cover for a given mass bin is highlighted by the colored hatches in panels a-c. Panel d overplots the r e — z distribution of the three mass bins 
(combining both barred and unbarred galaxies) and their lines of maximum SDSS fiber physical radius in their respective colors; the line of the minimum SDSS 
fiber physical radius is shown as the black horizontal line. The maximum SDSS fiber physical radius is « 50 — 70% of the median r e of each mass bin. 


range that is covered by the SDSS fiber at a given mass bin, 
which, as shown in Fig. [TJ corresponds to different maximum 
redshifts. Panel d overplots all the r e — z distributions of the 
three mass bins (combining both barred and unbarred galax¬ 
ies) and their lines of maximum SDSS fiber physical radius in 
their respective colors. 

Fig. |4] shows that for each mass bin, the maximum SDSS 
fiber physical radius is ss 50-70% of the median r e , mean¬ 
ing that the gradients in Fig. [3] extend out to « 50-70% the 
median galaxy semi-major half-light radius. Thus even at the 
highest redshift bin, our gradients are still dominated by the 
central regions of our galaxies, indicating that our gradient 
analysis is not sensitive to the outskirts of our galaxies. 

4.2.2. SDSS fiber size vs. typical bar length 

A similarly important question is “how much of our bars is 
covered by the SDSS fiber?” To illustrate the varying extent of 
the SDSS fiber with respect to the bar structure as a function 
of redshift, we present SDSS /'-band images of typical barred 
galaxies for each redshift interval at a given stellar mass bin 
with the SDSS fiber overplotted in Fig. [5) The physical size of 
the 3" SDSS fiber is shown at the bottom right of each image. 
For each mass bin, we reproduce the [Fe/H] gradients from 
Fig. Q] tor the reader’s convenience. 

Clearly, the SDSS fiber covers more of the bar with increas¬ 
ing redshift. This effect is not likely due to an evolution of 
the bar length: matching our barred galaxies from the gradi¬ 
ent analysis to the bar length catalog of iHovle et al.l (1201 Ill'll 
shows that there is no strong bar length evolution, indicating 

18 « 14% of our bars in the gradient analysis have bar length measure¬ 
ments from lHovle et alJ 1201 111 . This small overlap is because most galaxies 


that the SDSS fiber is truly covering a larger fraction of these 
bars with redshift. 

To answer the question posed at the beginning of this sec¬ 
tion more quantitatively, we need bar length measurements 
for all our barred galaxies. Unfortunately, as was mentioned 
above, we only have bar lengths for « 14% of our sample. 
Fortunately, this subsample is approxi matel y normally dis¬ 
tributed across the mass range used in 34.21 thus it is likely 
to be a representative sample. Assuming that this subsample 
is indeed representative of our entire bar sample, we find that 
most of our bars are 3-7 kpc (in semi-major axis), which 
is consistent with the findings of o ther works (lErwinl 120051 : 
iDurbala et al.1 l2008t iGadottJl20llh . Thus the SDSS fiber, 
which ranges from 0.6 kpc to 2.4 kpc in radius, covers « 8- 
80% of the length of the bars in our sample, i.e., the fiber 
samples a large range of the typical bar lengths in our sample. 

But it is clear that our gradients are limited to well within 
the bar structure, which is an important caveat to consider and 
is discussed in more detail in jj6] 

4.2.3. Stellar population gradients of barred and unbarred galaxies 

Given that our gradients do not extend beyond the bars, here 
we present the stellar population gradients of barred and un¬ 
barred quiescent disk galaxies in Fig. [3] The top three rows 
plot [Fe/H], [Mg/Fe], and [N/Fe] as a function of the physical 
radius of the SDSS fiber (in kpc) in the three bins of stellar 
mass listed above. The fourth row of Fig. [3] superimposes the 
gradients of all mass bins. At all radii, [Fe/H], [Mg/Fe], and 

in lHovle et alJ (2011ft have significant Ho emission, and moreover, they only 
considered pbar > 0.8 galaxies. 
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Figure 5. An illustration of the relative size of the SDSS fiber to the bar structure as a function of redshift. SDSS r-band images of typical barred galaxies at 
each redshift interval for a given stellar mass are shown at the top of their corresponding [Fe/H] gradient, which is reproduced from Fig. [3] Each image lists 
the galaxy’s SDSS DR7 object ID, stellar mass, debiased bar vote fraction, and redshift. The SDSS fiber is represented by the black circle at the center of each 
image with the correct relative scaling; the scale at the bottom right of each image gives the physical size of the 3” diameter of the SDSS fiber. The images are 
arranged from lowest redshift (i.e., smallest radius; left) to highest redshift (i.e., highest radius; right). Although the SDSS fiber covers more of the bar structure 
with redshift, the fiber is still within the extent of the bar at all redshifts. 


[N/Fe] increases with log M*, which is consistent with the 
bulge trends in Fig. [2] 

We stress that the stellar population parameters at a given 
fiber radius in Fig. |3]are not annular measurements, but rather, 
integrated measurements, i.e., we model the stellar popula¬ 
tion of the total area covered by the SDSS fiber, which in¬ 
creases with redshift. Despite the fact that our method of es¬ 
timating gradients is different from most other works (which 
use annular measurements), we still find results that are sim¬ 
ilar to those works. For example, our [Fe/H] gradients are 
negative (decreases with radius) while our [Mg/Fe] gradients 
are nearly-flat, both of which are consistent with the findings 
of past works dMehlert et all 120031: ISanchez-Blazquez et all 


2007 

:lRawle et al.! 2008:lKuntschner et alj 2010; Greene et al. 

2013 

). 


We do not consider the stellar age gradients because in addi¬ 
tion to probing larger physical distances within galaxies with 
redshift, we are also probing younger galaxies with redshift. 
Thus the resulting stellar age gradients can be due to both 
younger stellar populations at larger physical distances and 
passive evolution, a degeneracy that we do not attempt to dis¬ 
entangle. 

The main result of this section is evident upon examination 
of Fig. [3] there are no differences in the central stellar popula¬ 
tion gradients of barred and unbarred quiescent disk galaxies 
that exceed the level of systematic uncertainty. 
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5. COMPARISON TO PREVIOUS WORKS 


As summarized in the introduction, previous works 
in this topic often disagree—some show that bars 
affect their host galaxies’ stellar popukuions (e.g.. 


Moorthv & Holtzman 2006 

Perez & Sanchez-Blazquez 

2011; Sanchez-Blazquez et al. 

201 1: Coelho & Gadotti 2011; 

Williams et alj2012 

) while others show that bars do not (e.g.. 

Cacho et al. 2014; 

Sanchez-Blazquez et al. 2014j). While 


some of these differences can be attributed to the choice 
of stellar population synthesis model, IConrov et al.l (12014)) 
have shown that there is a broad agreement between their 
model (the one used in this wor k) and several other popula r 
models, including EZ_Ages (iGraves & Schiavoni 12008), 
which uses Lick indices as opposed to full-spectrum fitting. 
Thus assuming that different models produce similar results, 
we will restrict our comparison to data quantity, data quality, 
and sample selection. 

Comparing our data quantity to previous works shows that 
there are only a couple of works that are comparable— 
iCoelho & Gadottil (1201 1 | ) have 251 barred galaxies and 324 
unbarred galaxies, and ICacho et al.l (1201 4ft have 414 barred 
galaxies and 1180 unbarred galaxies; both works primarily 
contained star-forming barred galaxies. Our work contains 
over a 1,000 quiescent barred galaxies and over 5,000 quies¬ 
cent unbarred galaxies. 

Comparing our data quality to previous works shows that 
our stacked spectra have the highest S/N (> 100 A -1 ). Previ¬ 
ous works have spectra with S/N= 10-50 A -1 , which produces 
less accurate and less precise stellar population parameters. 
Of course, in order to achieve such high S/N, we have stacked 
many spectra, which has not been done in any previous work 
on this topic. Thus our work is measuring the average effects 
of bars in quiescent disk galaxies. 

Finally, comparing our sample selection to previous works 
shows that our sample is unique. While we only select galax¬ 
ies with SDSS spectra that are bereft of two of the strongest 
optical emission lines—Ha and [O ll]A3727—all previous 
works selected galaxies with a range of star formation states. 
In fa ct, m ost works require the presence of these lines (e.g., 
iCacho et all 120141) . The main benefit of our selection cri¬ 
teria is that our stacked spectra have minimal emission line 
contamination, resulting in highly accurate and precise stellar 
population modeling, which was one of the main goals of this 
work. This sample selection, however, is an important con¬ 
sideration in our interpretation, which we discuss in ^6] (see 
also !2. Ik 

Based on these three comparisons, it is clear that our work 
is unlike any of the previous works on this topic. Interest¬ 
ingly though, even with our differences, we find similar re¬ 
sults to the latest work on this topic. Namely, like lCacho et alJ 
d2014l) . we find no strong differences in the stellar populations 
of barred and unbarred galaxies, indicating that bars do not 
significantly affect the stellar populations of their host galax¬ 
ies, regardless of their star formation state. 

6. DISCUSSION 

6.1. Comparison to predictions of bar-driven secular 
evolution 

The results presented in this work indicate that bars do 
not affect—at least not strongly enough for us to signifi¬ 
cantly detect—the central stellar populations of quiescent disk 
galaxies. In the following discussion, we assume that bars are 
long-lived structures, which is consistent with the latest sim¬ 


ulations of bar formation and evolution (Athanassoul a et all 
l2013al) . Moreover, given that our sample comprises solely of 
quiescent galaxies, we will restrict our comparison to simu¬ 
lations without star formation. According to these types of 
simulations from [Friedli et al.l ( 1994 ). barred galaxies should 
have a lower central stellar [Fe/H] (by about 0.20 dex) com¬ 
pared to unbarred galaxies because the bar-driven inward 
transportation of stiy's^reduces the initial negative [Fe/H] gra¬ 
dient (see also lDi Matteo et ahll2013[ who find a reduction of 
about 0.10 dex in the central [Fe/H] of barred galaxies, albeit 
with unrealistic initial conditions). But as shown in Fig. [2] 
the bulge [Fe/H] of barred galaxies do not differ from that of 
unbarred galaxies at the few hundredths of a dex level. 

Chemical simulations of bar-driven secular evolution with¬ 
out star formation also predict that the slopes of the [Fe/H] 
gradients of bar red galaxies are flatter than those of un - 
barred galaxies dFriedli et al.l 119941 : iDi Matteo et alJ 1201 3ft . 
But Fig. 0 shows no obvious differences in the [Fe/H] slopes 
of barred and unbarred galaxies. It is important to note, how¬ 
ever, that our gradients do not extend beyond the bar, and are 
in fact, well within the bars of our sample. This is impor¬ 
tant because most simulations have shown that the flatten¬ 
ing of abundance gradie nts is most dramatic in the outskirt s 
of galaxies (e.g.. iFriedli et al 11 19941 IDi Matteo et al.ll2013ft . 
Nevertheless, simulations predict that the [Fe/H] gradients 
near the centers are _still affected by the presence of a bar— 
IDi Matteo et alJ (12013ft predicts a flattening by ~0.01 dex/kpc 
in the inner regions of their barred galaxy simulation—an ef¬ 
fect that we do not observe, albeit this effect may be too subtle 
for us to detect. 

The above mentioned simulations, however, do not include 
the effects of initial star formation and feedback, which are 
necessary for the accurate calculation of metallicity gradi¬ 
ents. Additionally, quiescent galaxies are not produced in any 
simulations, i.e., galaxies are either always star-forming or 


of bars is around z ~ 1 (e.g., Sheth et al. 

2008 

; Kraljic et ah! 

20121 Melvin et al. 2014]; Simmons et al. 

2014 

), and assum- 


ing that bars are long-lived (e.g.. lAthanassoula et alJl2013ah 

the simulations should use initial conditions based on what we 
know for galaxies of that era. No observational data are avail¬ 
able for metallicity gradients at such redshifts, apart from the 
work of lOuevrel et al.l d2012l) . which indicates that some of 
the gas-phase metallicity gradients of galaxies around z ~ 1.2 
could be positive even though most of them are negative. Fur¬ 
thermore, there are no corresponding measurements for the 
stellar metallicity gradients. Finally, these simulations men¬ 
tioned above assume smal l bulges, which may op timize the 
tangential force from bars (lLaurikainen et al.ll2007l) . Thus the 
predictions from these simulations may be overestimating the 
influence of bars. In view of this situation, it does not seem 
possible to compare adequately theoretical predictions with 
our results concerning the differential effect between barred 
and unbarred quiescent galaxies at z ~ 0. 

7. CONCLUSION 


In this work, we use the latest stellar population synthesis 
model of IConrov & van Dokkumld2012al) to estimate the stel¬ 
lar populations of barred and unbarred quiescent disk galax¬ 
ies with very high S/N (> 100 A -1 ) stacks drawn from the 
Sloan Digital Sky Survey. Taking advantage of the variable 
physical size of the SDSS fiber as a function of redshift, we 
model the stellar populations at different parts within galaxies. 
We first focus on the bulges of galaxies, which have a typical 
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half-light radius of 1 kpc (iFisher & DrorvlUolOh . and approx¬ 
imates the physical size of the SDSS fiber at 0.02 < z < 0.04. 
We then estimate the stellar population gradients by assuming 
that galaxies at a fixed stellar mass do not significantly evolve 
in the redshift range of 0.020 < z < 0.085, which spans a: 1 
Gyr. Our resulting stellar population gradients extend out to 
50%-70% of the median semi-major half-light radius of our 
samples, which is well within the extent of the bars in our 
sample. 

We find that there are no significant differences in the stel¬ 
lar populations of the bulges or the gradients of barred vs. 
unbarred quiescent disk galaxies, suggesting that bars are not 
a significant influence to the chemical evolution of quiescent 
disk galaxies. 
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APPENDIX 

UNBARRED QUIESCENT DISK SELECTION 

Face-on quiescent disk galaxies, especially those without structures such as bars, are difficult to distinguish from pure elliptical 
galaxies. In order to illustrate that our unbarred quiescent disk sample is not severely contaminated by pure elliptical galaxies, in 
this section of the appendix, we compare a sample of spheroid-dominate d gal axies to our unbarred quiescent disk s amp le. 

To select spheroid-dominated galaxies from our quiescent sample (see 32. lb we use the criteria outlined bv lCheng et all (1201 lb . 
In their work, they offer a way of selecting early-type, bulge-dominated galaxies without s harp e dges, which we consider a sat¬ 
isfactory definition of spheroid-dominated galaxies. The selection method of lCheng et al.l ( 201P depends on GIMD bulge+disk 
decompositions and standard SDSS photometric parameters. In particular, the criteria are: 

1. B/T >0.5 

2. s2 > 0.08 

3. b/a > 0.65 

4. 0 2.9 

with B/T representing the bulge-to-total ratio, s2 representing the smoothness parameter, b/a representing the r band axis 
ratio, and C representing the r band concentration. B/T and s2 are from the n = 4 + n = 1 GIM2D catalog (ISimard et al.l(2011 ) 
while b/a and C are from SDSS DR7. 

With this spheroid sample, we compare its distributions of global Sersic index (ri), central surface stellar mass density (E j kpc ; 
ICheung et al.l2012h . central velocity dispersion (a) , and bulge-to-total ratio (B/T) to that of our unbarred quiescent disk sample at 
the eight mass bins used in our bulge analysis (see 34. Il l in Figs. [6]and|7] We also include the distributions of our barred quiescent 
disk sample for a further comparison. 

These distributions show that the spheroids have, on average, a higher n, Ej kpc , er, and B/T than both the barred and unbarred 
quiescent disks for the four least massive bins, indicating that our unbarred quiescent disks are not dominated by spheroid- 
dominated galaxies. 

For the more massive bins, however, their distributions become more similar, indicating that the structural properties of these 
classes of galaxies are almost indistinguishable at high masses. This similarity though, does not mean that our unbarred quiescent 
disk sample is severely contained by ellipticals. Indeed, the structural distributions of our barred quiescent disk sample are also 
similar to that of the spheroid-dominated sample at these massive bins. And since bars are pure disk phenomena (lAthanassoula! 
1201 3hl) . we conclude that the structural parameters we considered simply cannot distinguish spheroid-dominated galaxies from 
quiescent disk galaxies at high stellar masses. Therefore, given the ineffectiveness of these structural parameters to distinguish 
face-on quiescent disks from spheroid-dominated galaxies at high masses, we rely on the visual morphologies of Galaxy Zoo. 

Perhaps the most reassuring result of this section, however, is that the structural distributions of our barred and unbarred 
quiescent face-on disks are similar, indicating that we are selecting similar types of galaxies. 
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Figure 6. Distributions of Sersic index ( n ), central surface stellar mass density (S* kpc ), central velocity dispersion (cr), and bulge-to-total ratio (B/T) for the 
barred quiescent disks, unbarred quiescent disks, and spheroid-dominated galaxies for the four least massive bins in the bulge analysis (see ED- 
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Figure 7. Same as Fig.[7]except for the four most massive bins in the bulge analysis. 
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Figure 8. Comparison of barred and unbarred stacked spectra for the four least massive bins in the bulge analysis (see E). The ratio of the barred to unbarred 
stacks are shown below the stacked spectra, with the grey region representing its error; several Lick indices highlighted. 

COMPARISON OF BARRED VS. UNBARRED STACKS 

In this section we directly compare the stacked spectra of the barred and unbarre d sam ples. Figs. [8]and[9]plot the barred stacks 
over the unbarred stacks for each stellar mass bin analyzed in the bulge analysis ( 34. 11 1. The bottom panels show the ratio of the 
barred to unbarred stacks, highlighting some relevant Lick indices. 

Clearly, there are very small differences in these stacks. To summarize these differences, we plot several Lick indices vs. stellar 
mass in Fig. [TO] H/3, a tracer of the stellar_age, [MgFe]’, a tracer for metallicity dThomas et al .1120031) , and Mg b and Mg b/<Fe>, 
a tracer for the o-abundance dThomas et akll2003i) all show no significant differences between barred and unbarred galaxies, 
mirroring our main result in 34.11 

COMPARISON OF MODEL FITS TO INPUT SPECTRA 

In this section of the appendix, we illustrate the quality of the model fits by sh owing the input stacked spectra and the model 
spectra in Figs. [TU and Q2]for four mass bins from the bulge analysis (see 34.1k the barred stacks are on the left and unbarred 
stacks are on the right. At the bottom of each spectra, we show the residuals of the model and input spectra. The first feature to 
note is that the residuals are small, about one percent, indicating that the fits are excellent. Secondly, the residuals for the barred 
and unbarred model spectra are similar, indicating that the quality of the fits of the two samples are comparable. 
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